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Complementary  metal-oxide  semiconductor  (CMOS)-hybrid  arrays  have  become  competitive  optical  de¬ 
tectors  for  use  in  ground-  and  space-based  astronomy.  Interpixel  capacitance  (IPC)  is  one  source  of  error 
that  appears  in  most  CMOS  arrays.  In  this  paper,  we  use  a  single-pixel-reset  method  to  model  IPC.  We 
combine  this  IPC  model  with  a  model  for  charge  diffusion  to  estimate  the  total  crosstalk  on  H4RG-10 
arrays.  Finally,  we  compare  our  model  results  to  55Fe  data  obtained  using  an  astrometric  camera  built  to 
test  the  H4RG-10  B0  generation  detectors.  ©  2012  Optical  Society  of  America 
OCIS  codes:  040.0040,  040.1240,  040.3060,  040.5160,  040.6040,  040.7480. 


1.  Introduction 

Complementary  metal-oxide  semiconductor  (CMOS) 
sensors  have  become  a  competitive  astronomical 
ground-  and  space-based  detector  solution.  CMOS 
sensors  have  a  flexible  readout  structure  that  allows 
a  single  pixel  or  group  of  pixels  on  the  array  to  be 
read  out  or  reset  at  any  time  without  disturbing  or 
reading  out  the  rest  of  the  array  This  random-access, 
nondestructive  read  capability  is  ideal  for  dynamic 
range-driven  astronomical  applications,  since  it 
permits  bright  and  faint  objects  to  be  observed 
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simultaneously  using  a  single  detector.  In  addition 
to  the  flexible  readout,  CMOS  sensors  are  naturally 
less  sensitive  to  radiation  than  more  traditional  de¬ 
tectors  such  as  charge  coupled  devices  (CCDs),  since 
damage  to  one  pixel  in  the  array  does  not  adversely 
affect  subsequent  pixels  in  a  row  or  column  of  the  ar¬ 
ray.  This  inherent  radiation  hardness  is  particularly 
appealing  for  space-based  applications. 

While  the  readout  capabilities  of  CMOS  sensors 
are  ideal  for  a  variety  of  observing  strategies,  the  fill 
factor  of  each  pixel  is  significantly  reduced  because 
the  readout  circuitry  is  implanted  directly  on  the 
photodetector  material.  Janesick  et  al.  [1]  describe 
how  the  CMOS-hybrid  focal  planes  have  been  devel¬ 
oped  to  address  this  issue.  A  CMOS-hybrid  sensor  is 
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a  CMOS  device  or  readout  integrated  circuit  (ROIC), 
mated  with  a  layer  of  photodetector  material.  The 
two  layers  are  typically  joined  together  using  indium 
bump  bonds.  The  resultant  hybrid  sensor  chip  as¬ 
sembly  (SCA)  is  back-illuminated,  and  combines  the 
flexible  readout  of  the  CMOS  with  a  CCD-like  fill  fac¬ 
tor  of  100%.  The  United  States  Naval  Observatory 
(USNO)  has  used  CMOS-hybrid  detectors  to  take 
advantage  of  this  performance  and  flexibility,  as  part 
of  its  development  of  very  large  format  focal  plane 
technologies,  in  support  of  the  Joint  Milli-Arcsecond 
Pathfinder  Survey  (JMAPS)  astrometry  mission. 

The  USNO  has  been  testing  large  format  Teledyne 
Imaging  Sensors  (TIS)  H4RG-10  Hybrid  Visible  Sili¬ 
con  Imager  (HyViSI)  SCAs  since  the  development 
of  the  first  generation-Al  detector  in  2006,  as  de¬ 
scribed  by  Dorland  et  al.  in  2007  [2].  In  2008  USNO 
supported  development  of  a  second  generation-A2 
H4RG-10  with  significantly  lower  dark  current  than 
the  Al  predecessor,  again  described  by  Dorland  et  al. 
in  2009  [3]  A  third  generation  detector,  the  H4RG-10 
B0  detector,  was  fabricated  in  2010  with  JMAPS 
support  to  address  yield  and  pixel  operability  re¬ 
quirements.  This  third  generation  detector  has  very 
low  dark  current  and  noise  properties  comparable  to 
many  CCDs.  The  quantum  efficiency  is  in  excess  of 
80%  across  multiple  wavelengths,  and  the  nonlinear¬ 
ity  is  <1%.  Additionally,  with  JMAPS  risk  reduction 
support,  Hubbs  et  al.  [4]  describe  how  TIS  was  able  to 
increase  pixel  operability  to  better  than  99.9%.  Here, 
operability  specifically  refers  to  the  percent  of  pixels 
that  are  fully  connected  without  shorting.  (For  a 
summary  of  the  development  of  the  H4RG-10  SCA, 
see  also  [5].) 

While  these  flexible  and  low-noise  CMOS-hybrid 
arrays  are  excellent  for  most  astronomical  applica¬ 
tions,  the  H4RG-10  shows  higher  levels  of  crosstalk, 
or  more  specifically,  the  interpixel  capacitance  (IPC) 
component  of  crosstalk,  than  CCDs.  We  also  note  that 
some  CMOS  pixel  circuit  designs  do  not  exhibit  high 
levels  of  IPC.  However,  the  chosen  design  of  the 
H4RG-10  has  advantages  that  are  difficult  to  achieve 
with  these  other  designs,  namely  low  read  noise,  low 
dark  current,  and  low  power  consumption.  Indeed, 
most  low  background  applications  prefer  the  source- 
follower  approach  despite  the  IPC  problem  (e.g., 
JMAPS  H4RG-10  B0),  while  higher  background  appli¬ 
cations  typically  use  one  of  the  other  design  options. 

As  discussed  in  detail  below,  crosstalk  can  be  pro¬ 
blematic  for  applications  like  astrometry,  because  it 
has  the  effect  of  blurring  the  point  spread  function 
(PSF)  of  the  photon  source,  resulting  in  lower  effective 
signal  to  noise  (S/N)  in  each  pixel  and  high  centroiding 
errors  (see  also  [6-8]).  For  this  reason  JMAPS  has 
supported  development  of  a  fourth  generation  detec¬ 
tor,  the  H4RG-10  Bl,  designed  in  part  to  reduce  IPC 
to  values  that  are  negligible  for  astrometry. 

In  this  paper,  we  discuss  the  data  analysis  and 
modeling  that  have  been  used  to  understand  cross¬ 
talk  for  astronomical  detectors.  In  Section  2,  we  de¬ 
scribe  the  models  for  charge  diffusion  and  IPC  that 


were  combined  to  create  a  55 Fe  model  based  on  the 
H4RG-10  B0  detector  (third  generation).  In  Section  3, 
we  describe  the  Jet  Propulsion  Laboratory  (JPL) 
camera  used  to  collect  the  data.  In  Section  4,  we  dis¬ 
cuss  the  analysis  methodology  for  single-pixel-reset 
(SPR)  data  and  compare  the  results  with  the  mod¬ 
eled  data.  We  also  discuss  the  analysis  methodology 
used  for  55  Fe  data  obtained  with  the  astrometric 
camera,  and  again  compare  the  analysis  and  model¬ 
ing  results.  Finally,  in  Section  5,  we  summarize  our 
findings.  (In  a  subsequent  paper,  we  show  how  cross¬ 
talk  affects  astrometry  and  photometry  using  de¬ 
tailed  detector  and  optics  models  for  a  realistic 
astrometric  telescope  (see  [6]  for  details.)) 

2.  Crosstalk  Models 

A.  Introduction  to  Charge  Diffusion  and  IPC 

Crosstalk  between  pixels  is  caused  by  two  indepen¬ 
dent  phenomena,  charge  diffusion  and  IPC.  Charge 
diffusion  is  the  lateral  movement  (i.e.,  pixel-to-pixel) 
of  charge  between  the  points  of  charge  production 
and  charge  collection  in  the  bulk  substrate  of  the  de¬ 
tector.  Charge  diffusion  occurs  in  all  photosensitive 
material,  including  the  silicon  substrate  of  all  CCDs. 
IPC  is  the  capacitance  that  arises  between  adjacent 
detector  pixels  in  the  source-follower  CMOS  design, 
and  leads  to  coupling  of  signal  between  those  pixels 
via  displacement  currents  flowing  from  the  collec¬ 
tion  node. 

For  optical  detectors  with  photosensitive  silicon 
material,  crosstalk  can  be  measured  by  exposing  the 
detector  to  a  radioactive  source  with  a  known  charge 
production  rate  in  bulk  silicon,  and  then  measuring 
the  voltage  in  pixels  surrounding  a  central  hit.  55  Fe 
is  a  soft  X-ray  source,  quickly  decaying  (half-life:  2.7 
years)  to  Mn  when  a  K  shell  electron  is  absorbed  into 
the  nucleus.  An  electron  from  either  the  L  or  M  shells 
drops  to  fill  the  hole  created  in  the  K  shell.  This  drop 
to  a  lower  energy  level  causes  the  emission  of  either  a 
Ka  (5.9  KeV)  or  Kp  (6.5  KeV)  X-ray.  The  absorption  of 
a  Ka  photon  in  the  bulk  silicon  of  the  detector  pro¬ 
duces,  on  average,  1620e“,  while  Kp  photon  absorp¬ 
tion  produces  an  average  of  1778c“. 

In  55 Fe  testing,  the  detector,  under  controlled  en¬ 
vironmental  conditions,  is  briefly  exposed  to  the  ra¬ 
dioactive  source,  and  the  resulting  charge  is  collected 
and  read  out  for  each  exposure.  When  analyzed, 
photon  hits  approximately  normal  to  the  detector 
surface  and  centered  on  a  single  pixel  are  averaged 
and  normalized  to  create  a  kernel  representative  of 
pixel  crosstalk.  55Fe  testing  and  analysis  is  described 
in  more  detail  in  Section  5. 

In  addition  to  55Fe  testing,  SPR  testing  allows  for 
the  direct  characterization  of  IPC  alone,  requiring  no 
illumination  source  [8] .  In  SPR,  after  setting  all  pix¬ 
els  in  the  SCA  to  a  single  voltage  and  making  an 
initial  readout,  a  well-spaced  grid  of  single  pixels 
is  reset  to  a  second  voltage  level.  A  subsequent  read¬ 
out  of  the  pixels  will  reveal  any  IPC  as  signal  in  pix¬ 
els  adjacent  to  the  reset  pixels.  As  with  55Fe  testing, 
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multiple  SPR  test  results  are  averaged  (discarding 
any  bad  pixels,  edge-affected  pixels,  etc.)  to  create 
a  representative  detector  IPC  kernel.  SPR  testing 
and  data  analysis  methodology  is  described  in  more 
detail  in  Section  4. 

The  IPC  kernel  obtained  through  SPR  testing, 
when  combined  with  charge  diffusion  modeling,  can 
be  used  for  55Fe  test  verification  and  prediction.  This 
resultant  55 Fe  model  is  based  on  the  convolution  of 
a  charge  diffusion  kernel,  resulting  from  an  incident 
Ka  photon,  with  an  IPC  kernel  built  on  SPR-measured 
IPC.  This  new,  modeled  kernel  can  be  used  to  approx¬ 
imate  the  average  pixel  crosstalk  expected  in  55Fe 
testing  for  a  given  detector.  The  subsections  below 
describe  the  modeling  of  IPC  and  charge  diffusion. 

B.  IPC  Model 

A  simple  model  of  the  IPC  expected  within  a  detector 
array  can  be  developed  based  on  the  approach 
described  by  Moore,  et  al.  in  2004  [9].  The  detector 
array  of  photodiodes  is  modeled  as  an  array  of  capa¬ 
citors,  as  shown  in  Fig.  1,  each  identical  and  with  a 
capacitance  that  is  unchanging  with  voltage. 

The  detector  is  constructed  in  a  way  that  allows  an 
electrical  field  to  exist  between  neighboring  collec¬ 
tion  nodes,  essentially  creating  small  coupling  capa¬ 
citors  between  the  nodes.  Charge  entering  a  single 
nodal  capacitor  Qtotal  causes  a  voltage  change  in  that 
node,  and  through  the  coupling  capacitors,  causes 
voltage  changes  in  neighboring  nodes.  Qtotai  is  the 
sum  of  all  the  apparent  charge  seen  both  in  the  vol¬ 
tage  of  the  original  node  and  the  voltages  of  n  neigh¬ 
boring  nodes. 


Qtotai 

^node 


(1) 


In  Moore’s  approach,  the  impulse  response  of  each 
node  h(n)  is  a  ratio  of  the  charge  that  appears  elec¬ 
trically  in  a  node  Qn ,  to  the  photocurrent  that  en¬ 
tered  the  original  node  Qtotai? 


h(n)  = 


Qn 

Qtotai 


(2) 


Fig.  1.  (Color  online)  Interpixel  capacitance  (Cip). 
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Fig.  2.  IPC  model  for  3  x  3  pixel  array. 


Using  this  approach,  and  assuming  two- 
dimensional  symmetry  in  nearest  neighboring  pixels 
and  diagonal  neighboring  pixels,  a  simple  model  of 
IPC  can  be  constructed.  Charge  appearing  electri¬ 
cally  in  surrounding  pixels  is  defined  in  terms  of  a 
single  variable  a,  defined  as  the  percent  of  total 
charge  seen  in  any  of  the  four  nearest  neighbor  pix¬ 
els.  Symmetry  in  nearest  neighbors  and  diagonal 
neighbors  is  assumed  in  this  model.  An  example  is 
shown  in  Fig.  2  for  a  3  x  3  pixel  array. 

The  sum  of  the  charge  in  all  pixels  in  the  array  = 
1.0.  The  above  model  kernel,  with  a  =  0.075,  is  an 
IPC  model  for  a  general  CMOS  detector.  The  model 
can  be  extended  to  larger  kernel  models  to  account 
for  detectors  with  IPC  that  is  more  broadly  spread. 

C.  Charge  Diffusion 

Lateral  charge  diffusion  occurs  while  charge  moves 
between  the  point  of  generation  and  the  point  of  col¬ 
lection,  in  the  detector  substrate.  The  process  begins 
with  the  absorption  of  a  Ka  photon,  which,  for  55 Fe, 
produces  a  cloud  of  1620  charge  pairs  within  the  sub¬ 
strate,  as  illustrated  in  Fig.  3. 

The  photon  attenuation  length  la  for  a  Ka  (5.9  KeV) 
photon  is 


a  4n  Im(72)  ’  1  J 

using  the  absorptive  term  Im(72)  of  the  complex 
refractive  index  of  silicon,  n, 


Detector  backside 


V=A0\ 


Si  substrate 


Charge  cloud 


z 

A 


y 


zD  =  1 00  um 


Fig.  3.  (Color  online)  Charge  diffusion  decreases  with  increasing 
absorption  depth  along  z. 
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n  =  1  -  8  -  i/3, 


(4) 


where  8  is  the  refractive  index  decrement  and  (3  is 
the  absorptive  index.  For  a  photon  with  E  =  5.9  KeV, 
~  6.  Ox  10-7  [10]  and  la  ~  28  / <m .  In  the  55Fe  model, 
charge  diffusion  is  calculated  for  thin  slices  of  the 
bulk  Si  substrate.  The  probability  of  photon  absorp¬ 
tion  is  calculated  for  each  slice: 


E 


max 


qNd 

2esi 


When  V  is  much  larger  than  VD  [12,16], 


a(z)  ~ 


2z%kBT 
qV  ’ 


(ID 


(12) 


1  /  ~^1  ^2  \ 

Paz  =  T{e^-e^\,  (5) 

where  di  and  d2  are  positions  along  the  z-axis  (a  di¬ 
rection  normal  to  the  pixel  array  plane  ( xy  plane),  as 
defined  for  illustration  in  Fig.  3  above),  and  d1  is 
further  from  the  detector.  (In  this  model,  the  incident 
photon  is  always  normal  to  the  detector  surface.  This 
orientation  was  selected  to  mimic  the  symmetry  se¬ 
lection  criteria  for  good  hits  used  in  actual  55Fe  test 
data  analysis,  as  described  in  Subsection  4.D.) 

We  now  consider  charge  diffusion  for  a  detector 
with  a  fully  depleted  substrate.  Charge  diffusion,  as 
a  function  of  depth  and  weighted  for  absorption  prob¬ 
ability,  is  well  represented  by  the  two-dimensional 
Gaussian 


f  CD  ~ 


Pa  z 

d2na(z) 


g  2 a(zY 


(6) 


with  <r(z)  the  root  mean  square  (RMS)  standard  de¬ 
viation  of  charge  spreading  in  jc  andy.  This  spreading 
is  described  in  terms  of  the  diffusion  constant  Dp  and 
the  transit  time  ttY  for  charges  (holes,  for  the  detec¬ 
tors  under  consideration)  to  move  from  the  point  of 
charge  pair  generation  to  the  point  of  collection 


g(z)  =  j2Dpt,.  (7) 

The  transit  time  tir  can  be  found  using  the  hole 
drift  velocity  ydrift  as  described  by  [11,12] 


^  drift  ^  flE(z) 


(8) 


and  integrating  over  the  entire  depletion  depth, 


^tr 


€Si 

qNd 


(9) 


where  e^i  is  the  dielectric  constant  for  silicon  (11. 9x 
8.854  x  10-12  C/V),  q  is  the  fundamental  charge  unit 
(C),  and  Nd  is  the  doping  constant  for  the  bulk  silicon 
substrate.  (We  use  a  doping  density  of  1  x  1018  m-3, 
which  is  a  typical  doping  density  for  silicon  and  re¬ 
presents  the  average  of  values  presented  in  [13- 
15].)  Emax  includes  the  effects  of  both  the  substrate 
bias  voltage  V  and  the  depletion  voltage  VD  [16] , 


v  _QNd72 

Vd  ~  2iT°' 


(10) 


where  the  Einstein  relation,  DP//i  =  kBT/q ,  has  been 
applied  to  Eq.  (7),  above,  and  where  T  is  the  tempera¬ 
ture  of  the  detector  ( K )  and  kB  is  Boltzmann’s  con¬ 
stant  ( J/K ).  For  absorption  of  the  photon  at  a 
specific  z  within  the  substrate, 


o(z)  = 


2kBTeSl 
.  q2Nd 


In 


(  ffmax  V 

-2z-zd))1 


1/2 


(13) 


In  the  model,  the  charge  diffusion  function  fCD  is 
calculated  for  each  incremental  slice  through  z,  using 
the  above  equation  for  g(z),  with  T  =  193  K,  V  =  40  V, 
and  zD  =  100  //m. 

The  charge  cloud  produced  by  absorption  of  a  Ka 
photon  is  resolved  at  a  higher-than-pixel  resolution, 
and  diffusion  to  each  neighboring  pixel  is  determined 
by  the  position  of  the  cloud  center  projected  some¬ 
where  on  the  surface  of  the  central  pixel.  Because 
of  this,  the  charge  diffusion  is  calculated  at  a  subpixel 
resolution,  with  an  absorption  probability  weighting 
for  each  slice.  The  high-resolution  kernel  is  then 
summed  and  rebinned  to  a  lower  resolution  kernel, 
ready  for  convolution  with  the  IPC  kernel.  Examples 
of  the  charge  diffusion  model  to  both  subpixel  and 
pixel  resolutions  are  shown  in  Fig.  4. 

3.  JMAPS  Test  Camera 

The  data  described  here  were  taken  at  the  JPL  using 
a  test  camera  specifically  designed  for  ground  based 
astrometric  testing  of  H4RG-10  detectors  for  JMAPS. 
The  camera  functions  with  32  outputs  controlled 
by  a  TIS  noncryogenic  SIDECAR  (system  image, 
digitizing,  enhancing,  controlling,  and  retrieving)  ap¬ 
plication  specific  integrated  circuit  (ASIC).  The  mea¬ 
surements  were  taken  at  193  K,  and  the  detector 
substrate  voltage  was  set  to  10,  20,  30,  or  40  V,  de¬ 
pending  on  the  measurement.  The  camera  voltage 
settings  were  the  same  for  IPC  SPR  and  55 Fe  data 
acquisition.  Table  1  lists  the  primary  voltage  para¬ 
meters  used  for  these  measurements.  An  image  of 
this  camera  is  shown  in  Fig.  5. 

4.  SPR  Data 

A.  SPR  Data  Acquisition  and  Analysis 

The  JMAPS  astrometric  test  camera  is  equipped 
with  microcode  that  can  perform  the  SPR  technique 
described  by  Seshadri  [17]  and  Finger  [1_8] .  The  SPR 
function  enables  reset  of  single  pixels  in  the  array  to 
a  different  reset  voltage  than  the  rest  of  the  pixels  in 
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Fig.  4.  (Color  online)  Charge  diffusion  example  images,  with  Qpixei/Qtotai  for  each  pixel. 


the  array.  The  effective  result  is  charge  on  the  inte¬ 
gration  node  for  these  reset  pixels.  This  method  is 
useful  for  isolating  IPC,  since  the  charge  is  not  gen¬ 
erated  in  the  photodetector  material,  and  is  therefore 
not  susceptible  to  the  effects  of  charge  diffusion. 

IPC  maps  are  created  using  the  following  data  ac¬ 
quisition  method.  First,  the  entire  detector  is  reset 
to  a  given  reset  voltage  and  the  pixels  are  read  out 
to  generate  an  offset  frame.  Then,  a  grid  of  widely 
spaced  single  pixels  is  reset  to  a  new  voltage.  Follow¬ 
ing  this  reset,  all  pixels  are  again  read  out.  If  IPC  is 
present,  the  difference  of  these  two  images  will  show 
voltage  in  pixels  neighboring  the  reset  single  pixels. 
A  cropped  SPR  image  of  a  section  of  an  H4RG-10  BO 
detector  taken  with  the  JMAPS  camera  is  shown 
in  Fig.  6. 

For  this  analysis,  the  initial  reset  voltage  was  first 
set  to  300  mV  for  all  pixels.  The  SPR  grid  of  pixels 
was  then  reset  to  0,  100,  200,  300,  400,  500,  600  mV. 
This  resulted  in  seven  measurements  of  varying  sig¬ 
nal  levels  per  substrate  voltage  (10  V,  20  V,  30  V, 
40  V).  To  analyze  the  data,  the  following  procedure 
for  finding  good  hits  was  adopted. 

1.  First  the  image  was  offset  (bias)  corrected. 

2.  A  small  range  of  central  pixel  values  were 
identified  as  hits.  The  range  of  values  for  the  central 
pixel  was  identified  as  a  Gaussian  distribution  in  the 
histogram  of  the  image,  and  was  well  above  the  noise 
floor  of  the  image.  Defining  this  range  is  important  to 
ensure  that  hot  or  bad  pixels  are  excluded  from  the 


Table  1.  JMAPS  Camera  Voltage  Parameters 


Parameter 

Symbol 

Voltage  (V) 

Bias  voltage 

VSUB 

10-40 

Digital  positive  power  supply 

VDD 

3.12 

Analog  positive  power  supply 

VDDA 

3.12 

Drain  node  of  pixel  SF 

CELLDRAIN 

0.00 

Drain  node  of  output  SF 

DRA 

0.00 

Source  node  of  internal  current 
for  SF 

VBIASPOWER 

3.00 

Detector  substrate  voltage 

DSUB 

0.50 

Detector  reset  voltage 

VRESET 

0.32 

Bias  voltage  for  current  source 
of  SF 

VBIASGATE 

2.00 

IPC  estimate.  On  average,  each  image  contains  ap¬ 
proximately  64  000  good  hits  that  are  used  for  the 
analysis. 

3.  All  pixels  defined  as  hits  based  on  the  criteria 
above  were  background  corrected  using  a  local  back¬ 
ground  correction  method.  An  annulus  of  two  pixels 
outside  the  central  hit  was  extracted.  The  median  of 
the  pixels  in  this  annulus  was  subtracted  from  the 
values  of  the  inner  9x9  pixels.  An  image  illustrating 
the  annulus  and  kernel  is  shown  in  Fig.  7. 

4.  The  background-corrected  kernels  were  then 
stacked  and  averaged. 

5.  The  resultant  measured  IPC  kernel  was 
normalized,  so  that  the  charge  in  a  single  pixel 
was  represented  as  a  percent  of  the  total  charge  in 
the  kernel. 

For  this  detector,  the  noise  floor  was  reached  in 
both  IPC  and  55Fe  crosstalk  data  at  the  9x9  pixels 
boundary. 

B.  IPC  Data  and  Model 
1.  Data 

Figure  8  shows  the  fraction  of  signal  in  the  central 
pixel  as  a  function  of  the  SPR  voltage.  In  this  plot, 
the  saturation  limit  is  indicated  by  the  drop-off  in 
fractional  signal  at  400  mV.  However,  in  order  to  be 
certain  that  no  saturated  pixels  were  included  in  the 
final  averaged  kernel,  and  to  ensure  high  S/N,  the 


Fig.  5.  (Color  online)  JMAPS  camera  and  H4RH-10  B0  in  the 
dewar. 
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Fig.  6.  Image  of  pixels  that  have  been  reset  using  the  SPR  meth¬ 
od  for  estimating  IPC. 


200  mV  SPR  data  set  was  chosen  for  this  analysis  (as 
opposed  to  0,  100,  and  300  mV  SPR  data). 

Figure  9  shows  the  9x9  pixels  kernel  resulting 
from  the  data  taken  at  40  V  and  200  mV  SPR.  Despite 
the  fact  that  the  outer  9x9  annulus  in  the  kernel  be¬ 
low  appears  to  be  statistically  indistinguishable  for 
the  errors,  we  note  that  all  values  in  this  annulus  are 
systematically  above  zero,  indicating  that  the  noise 
floor  has  not  been  reached.  For  this  reason,  the 
9x9  annulus  is  included  in  the  kernel.  The  values 
listed  represent  the  percentage  of  the  total  signal 
for  an  average  kernel,  based  on  approximately 
64  000  events.  Measurement  errors  represent  the 
standard  deviation  across  all  SPR  events. 

2.  Model 

As  can  be  seen  in  the  IPC  data  kernel  in  Fig.  9,  IPC 
can  spread  significantly  beyond  the  3x3  grid 
depicted  in  Subsection  2.B.  In  fact,  for  each  detector 
tested  using  SPR,  the  IPC  distribution  extended  to  at 
least  a  9  x  9  grid.  Additionally,  each  detector  type 
(where  type  is  determined  by  changing  design 
and  processing)  was  found  to  have  a  unique  IPC 


0  2  4  6  a  to  u 


Fig.  7.  (Color  online)  Illustration  of  background  area  used  for  cor¬ 
rection.  The  shaded  region  of  2-pixel  width  is  the  annulus  used  for 
correction. 


SPR:  Fractional  Voltage  in  Center  Pixel  vs  V,^,  at  40V 


Fig.  8.  (Color  online)  Fractional  voltage  in  the  central  pixel  as  a 
function  of  the  reset  voltage  level  for  SPR  data. 


distribution  signature.  Based  on  these  findings, 
the  IPC  model  is  an  averaged  set  of  SPR  data  for  each 
detector  type,  and  given  in  terms  of  a  (the  average 
fractional  signal  in  nearest  neighboring  pixels,  per 
Subsection  2.B).  The  use  of  a  is  in  keeping  with 
the  standard  modeling  approach  for  IPC  and  pro¬ 
vides  a  variable  to  capture  the  standard  deviation 
of  the  measured  and  summed  SPR  kernels  in  the 
55Fe  model  (discussed  in  Subsection  4.E).  Variation 
of  a  is  also  needed  to  model  the  effects  of  varying  le¬ 
vels  of  optical  crosstalk  on  centroiding  precision. 

In  the  example  above  (Fig.  9),  a  =  (4.51+ 
4.61  +  4.61  +  4.51)/4.0  =  4.56.  Using  this  a,  an  IPC 
model  (Fig.  10)  was  developed  to  capture  the  signal¬ 
spreading  signature  of  the  B0  generation.  Each  cell  is 
a  factor  of  a ,  where  the  central  pixel,  C  =  (1/a)  J2  (°f 
the  grid  of  coefficients). 

C.  IPC  Asymmetry  Findings 

Given  the  ROIC  readout  asymmetry  in  these  detec¬ 
tors,  we  specifically  set  out  to  determine  if  a  corre¬ 
sponding  asymmetry  exists  in  the  IPC.  The  JMAPS 
test  camera  at  JPL  is  designed  to  read  out  the  full 
4096  x  4096  pixel  detector  using  32  equally  spaced 
outputs.  This  means  that  the  detector  is  virtually 
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Fig.  9.  (Color  online)  IPC  values  for  H4RG-10  B0  based  on  SPR 
method. 
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Fig.  10.  (Color  online)  IPC  model  based  on  SPR  data. 


divided  into  32  128  x  4096  pixel  columns.  The  col¬ 
umns  are  read  out  alternately  from  the  right  and  left, 
resulting  in  all  odd  columns  being  read  out  in  one 
direction  and  even  columns  in  the  other  direction. 
Figure  11  shows  this  left/right,  odd/even  dichotomy. 

We  conducted  a  column-by-column  comparison  of 
the  IPC  charge  distribution  for  BO  detectors.  To  per¬ 
form  this  analysis,  the  hits  from  all  even  and  all  odd 
columns  were  averaged  separately.  The  first  and  last 
columns  of  the  array  were  excluded  from  this  analy¬ 
sis,  since  they  suffer  high  dark  current  at  the  edges 
that  could  contaminate  the  results.  A  clear  horizon¬ 
tal  dichotomy  resulting  from  the  direction  of  the 
readout  is  visible  in  the  odd  versus  even  columns, 
as  shown  in  Fig.  12. 

We  performed  this  analysis  for  all  of  the  SPR  data 
taken  for  0-300  mV  Vreset  values  and  for  all  four  sub¬ 
strate  voltages.  Table  2  summarizes  our  findings.  In 
the  readout  direction  (horizontal  direction),  the  offset 
induced  by  the  readout  is  approximately  5%,  depend¬ 
ing  on  whether  the  column  is  odd  or  even.  In  addition, 
a  vertical  bias  is  evident  in  every  column  (both  even 
and  odd)  of  approximately  3.5%.  The  fifth  and  eighth 
columns  of  Table  2  show  the  residuals  resulting  when 
hits  from  both  odd  and  even  columns  are  averaged. 
We  find  the  residuals  to  be  low  in  the  horizontal/read 
direction  after  averaging  odd  and  even  columns,  but  a 
3.5%  residual  is  still  present  in  the  vertical  direction. 
These  findings  will  be  used  in  the  selection  criteria  of 
55 Fe  hits,  as  described  in  Subsection  4.D. 


IPC  Kernel  Model  (%  alpha) 


Fig.  11.  (Color  online)  Readout  “columns”  determined  by  alter¬ 
nating  readout  direction. 


“Odd”  Readout,  128  pixels  wide 


V 

“Even”  Readout,  128  pixels  wide 


Fig.  12.  (Color  online)  SPR-generated  IPC  kernels  from  even 
(left)  and  odd  (right)  columns. 


D.  55Fe  Data  Acquisition  and  Analysis 

55Fe  measurements  were  taken  at  the  same  voltage 
settings  as  the  SPR  measurements,  using  the 
JMAPS  camera.  For  these  measurements,  the  source 
was  exposed  for  up  to  12  seconds.  During  the  first  two 
seconds,  a  mechanical  shutter  shielded  the  detector 
from  the  source.  This  was  done  to  measure  the  offset 
from  the  zero  point,  due  to  reset,  without  contamina¬ 
tion  from  the  55Fe  hits.  In  the  remaining  10  seconds, 
the  shield  was  removed  and  the  detector  was  nondes- 
tructively  read  out  every  2  seconds,  resulting  in  five 
frames  of  a  progressively  increasing  population  of 
55Fe  hits. 

The  charge  diffusion  component  of  crosstalk  is 
heavily  dependent  on  the  depth  at  which  the  Ka  or 
Kp  photon  is  absorbed.  Two  things  must  be  taken 


Table  2.  Results  of  the  Odd/Even  Column  Systematic  Offset  Test 
13  cm 


Vsuba  (V) 

Vreset6  (V) 

xec 

x0d 

•Terror 

yj 

y08 

terror 

10.00 

0.30 

1.08 

0.93 

-0.23% 

0.94 

0.94 

6.01% 

10.00 

0.20 

1.07 

0.94 

-0.29% 

0.97 

0.97 

3.27% 

10.00 

0.10 

1.06 

0.94 

-0.28 

0.97 

0.97 

3.34% 

10.00 

0.00 

1.05 

0.96 

-0.39% 

0.96 

0.96 

3.52% 

20.00 

0.30 

1.08 

0.93 

-0.29% 

0.97 

0.97 

2.87% 

20.00 

0.20 

1.07 

0.94 

-0.28% 

0.97 

0.97 

3.16% 

20.00 

0.10 

1.06 

0.94 

-0.27% 

0.97 

0.97 

3.32% 

20.00 

0.00 

1.05 

0.96 

-0.42% 

0.97 

0.97 

3.39% 

30.00 

0.30 

1.08 

0.93 

-0.33% 

0.97 

0.97 

2.96% 

30.00 

0.20 

1.07 

0.94 

-0.26% 

0.97 

0.97 

3.27% 

30.00 

0.10 

1.06 

0.94 

-0.26% 

0.97 

0.97 

3.38% 

30.00 

0.00 

1.05 

0.95 

-0.34% 

0.97 

0.97 

3.43% 

40.00 

0.30 

1.08 

0.92 

-0.28% 

0.97 

0.97 

3.07% 

40.00 

0.20 

1.07 

0.93 

-0.21% 

0.97 

0.97 

3.32% 

40.00 

0.10 

1.06 

0.94 

-0.24% 

0.96 

0.97 

3.49% 

40.00 

0.00 

1.05 

0.95 

-0.30% 

0.96 

0.96 

3.53% 

Averages 

-0.29% 

3.46% 

“substrate  voltage. 

6SPR  reset  voltage. 

bright  Aieft  f°r  all  even  columns. 

bright  Aieft  for  all  odd  columns. 

e%  residual  error  from  hits  in  both  columns. 

(y upper / y lower  for  even  columns. 

Supper  /blower  for  odd  columns. 

h%  residual  error  from  hits  in  both  columns. 
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into  consideration  when  selecting  good  hits.  First, 
the  primary  goal  is  to  select  hits  that  are  as  perpen¬ 
dicular  to  the  detector  surface  plane  as  possible.  Sec¬ 
ond,  the  hits  must  be  selected  so  that  all  absorption 
depths  are  sampled.  We  define  an  accurate  measure¬ 
ment  of  55Fe  as  the  average  over  all  absorption 
depths.  To  ensure  this,  we  have  used  the  following 
method  to  define  normal  55 Fe  hits: 

1.  First  the  image  was  offset  corrected.  To  avoid 
any  noise  associated  with  settling,  the  frames  used 
for  this  analysis  were  the  fourth  and  fifth  frames. 
The  fourth  frame  was  taken  as  the  offset  frame  and 
subtracted  from  the  fifth  frame. 

2.  A  small  range  of  central  pixel  values  were 
identified  as  hits.  The  range  chosen  is  in  analog  di¬ 
gital  units  (ADUs)  and  strongly  depends  on  the  gain 
of  the  system.  Therefore,  for  this  initial  rough  esti¬ 
mate,  the  range  was  taken  directly  from  a  subsample 
of  hits  on  the  image  without  regard  to  the  symmetry 
of  the  hit.  The  sum  in  a  9  x  9  region  around  the  cen¬ 
tral  pixel  represents  the  majority  of  the  ADUs  result¬ 
ing  from  the  hit.  A  Gaussian  distribution  centered  on 
the  Ka  peak  results  when  these  sums  are  plotted  in  a 
histogram.  We  defined  the  median  of  a  Gaussian  fit 
to  this  distribution  divided  by  the  total  electrons  ex¬ 
pected  for  a  Ka  hit  (1620e“)  as  the  gain  of  the  system. 
In  our  case,  the  gain  was  1.19  +  /  -  0.04  ADUs  per 
electron.  On  average,  ~23, 000  hits  proceed  to  step  3 
as  potential  Ka  hits. 

3.  Double  hits,  or  hits  falling  inside  the  param¬ 
eter  of  the  targeted  hit  region,  can  skew  the  crosstalk 
results.  To  remove  double  hits,  all  of  the  potential  Ka 
hits  from  step  2  were  centered  and  stacked  in  the 
z  direction,  a  process  depicted  in  Fig.  13.  Kernels 
with  double  hits  were  removed  by  fitting  a  Gaussian 
distribution  along  the  z- axis  of  the  stack  for  each  pix¬ 
el  except  the  central  3  x  3  in  the  kernel  stack.  Hits  in 
the  stack  with  pixel  values  greater  than  the  mean 
+/  -  3cr  over  the  stack  were  considered  double  hits, 
and  the  hit  was  removed  from  the  analysis.  The  re¬ 
maining  hits  in  the  stack  are  considered  single  hits 
and  were  used  in  the  following  steps  of  the  analysis. 

4.  A  local  dark  correction  was  then  performed  on 
the  remaining  single  hits.  The  pixels  in  an  annulus  of 
2  pixels  outside  of  the  9x9  central  hit  are  defined  as 


Fig.  13.  (Color  online)  Double  hits  are  removed  by  fitting  a  Gaus¬ 
sian  to  stacked  pixels.  Pixels  with  crosses  were  excluded  from 
iteration. 


background  pixels  (see  Fig.  7).  The  median  of  the  pix¬ 
els  in  this  annulus  was  subtracted  from  the  values  of 
the  inner  9x9  pixels. 

5.  The  symmetry  criteria  are  perhaps  the  most 
challenging  criteria  to  incorporate  into  the  55 Fe  re¬ 
duction  method.  Noise  and  detector  effects  prevent 
even  a  perfectly  symmetric  hit  from  appearing  sym¬ 
metric  in  the  data.  This  leads  to  confusion  of  off- 
center  hits  with  perfectly  centered  hits  contaminated 
by  noise.  To  choose  perfectly  centered  hits,  we  use  the 
following  criteria: 

(a)  The  noise  of  the  system  can  be  measured  by  fit¬ 
ting  a  Gaussian  distribution  to  all  pixels  in  the  image 
except  those  affected  by  55Fe  hits.  The  sigma  of  this 
distribution  represents  the  1-sigma  noise  of  the  mea¬ 
surement.  Two  times  this  sigma  is  the  maximum  dif¬ 
ference  in  counts  that  the  left/right  and  upper/lower 
pixels  can  have,  due  to  the  noise.  For  all  of  our  mea¬ 
surements,  this  noise  was  ~15  ADUs  (1(7). 

(b)  As  discussed  in  Subsection  4.C,  the  B0  detec¬ 
tors  show  both  horizontal  and  vertical  systematic 
readout  asymmetry.  The  asymmetry  in  the  horizon¬ 
tal  direction  is  column-dependent,  while  the  asym¬ 
metry  in  the  vertical  direction  is  not.  We  used  the 
normalized  kernels  from  our  IPC  analysis  to  esti¬ 
mate  the  offset  in  ADUs  expected  to  result  from  the 
asymmetry  induced  by  the  readout.  In  most  cases, 
this  asymmetry  in  effect  contributed  ~5  -  10  ADUs 
to  the  total  noise. 

(c)  With  the  noise  and  asymmetry  tabulated,  we 
defined  centered  hits  as  those  hits  having  left/ 
right  and  upper/lower  differences  that  are  less  than 
2-sigma  plus  the  asymmetry  offset.  In  most  cases, 
this  resulted  in  a  total  permissible  variation  between 
left/right  and  upper/lower  pixels  of  ~35  -  40  ADUs. 

6.  Finally,  in  order  to  ensure  only  Ka  hits  are  in¬ 
cluded  in  the  final  crosstalk  kernel  (for  the  purposes 
of  comparing  the  data  to  the  model),  the  sum  over  all 
9x9  pixels  of  each  symmetric  hit  was  taken.  Using 
the  gain  from  step  2,  any  hit  with  a  sum  from  the 
9x9  that  had  values  less  than  or  equal  to  the  Kp  es¬ 
cape  peak  or  greater  than  or  equal  to  the  Kp  peak  [12] 
were  excluded  from  the  stack.  We  note  that  the  peaks 
of  the  Kp  lines  are  used,  rather  than  the  edges  of 
these  peaks  (based  on  Janesick  [12]),  in  choosing  our 
range.  We  do  this  because  the  noise  is  high  enough 
for  our  system  that  the  Kp  peaks  cannot  be  distin¬ 
guished  from  the  Ka  peak  when  plotted.  To  avoid  pre¬ 
ferentially  excluding  Ka  hits  that  overlap  with  the 
Kp  hits  due  to  noise,  we  established  this  definition. 
We  recommend  using  only  Ka  hits  for  those  systems 
with  noise  that  is  low  enough  to  distinguish  the  three 
peaks. 

7.  After  applying  these  criteria,  the  remaining 
hits  were  averaged  to  generate  a  9x9  kernel  (in 
the  case  of  the  40  V  data)  in  units  of  ADU.  The  cross¬ 
talk  kernel  was  defined  as  the  charge  in  a  single  pixel 
divided  by  the  total  charge  in  the  9x9  kernel. 

8.  This  test  was  repeated  nine  times,  resulting 
in  nine  kernels  for  all  nine  data  sets.  The  standard 
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deviation  across  all  nine  measurements  for  each 
pixel  represents  the  error  on  the  final  kernel  due 
to  statistical  differences  between  each  data  set. 


This  method  describes  the  reduction  strategy  for 
the  40  V  55Fe  data.  With  lower  substrate  voltage  set¬ 
tings,  we  find  that  a  kernel  larger  than  9  x  9  is 
needed  to  capture  all  of  the  escaped  charge  from 
the  central  pixel.  We  discuss  this  in  more  detail  in 
Subsection  4.F. 

E.  55Fe  40  V  Data  and  Model 

1.  Data 

Figure  14  shows  the  55Fe  kernel  resulting  from  the 
acquisition  and  analysis  strategy  described  in  Sub¬ 
section  4.D  for  the  H4RG-10  BO  detector.  Kernel  pixel 
values  are  the  percent  measured  signal  for  an  aver¬ 
age  of  nine  tests.  The  total  number  of  good  hits  used 
for  this  kernel  was  1094.  The  errors  are  also  shown 
for  each  pixel  value  and  represent  the  standard  de¬ 
viation  of  the  pixel  values  over  all  nine  data  sets.  We 
define  pixel  crosstalk  as  the  value  of  the  average  four 
nearest  neighbor  pixels  (shown  in  gray)  divided  by 
the  central  value  (shown  in  yellow).  For  the  H4RG- 
10  BO  detectors,  we  obtain  a  measured  crosstalk 
value  of  9.83%  for  this  55Fe  data  set. 


2.  Model 

To  generate  the  55Fe  model,  the  average  IPC  kernel 
(/’iPC,j/)>  as  described  above  in  Subsection  4.B,  is  used 
to  represent  the  IPC  spreading  typical  of  the  detec¬ 
tor.  The  uncertainty  associated  with  the  standard  de¬ 
viation  (<Tcp)  of  the  averaged  SPR  results  is  included 
in  the  55 Fe  model  through  convolution  of  the  charge 
diffusion  kernel  with  a  set  of  IPC  kernels  derived 
from  new  alphas  (crnew)-  These  anew  values  are  calcu¬ 
lated  from  a  set  of  central  pixel  values  (Cnew)  selected 
from  a  Gaussian  distribution  built  around  the  aver¬ 
aged  central  pixel  value  (Cave)  and  the  standard  de¬ 
viation  of  that  averaged  value,  <rcp: 


Cnew  =  ^cp[-2  In  (1  -  raj)]1/2  cos  (2 nrn2)  +  Cave,  (14) 


(i-cnew) 

J2ijf  IPC, ij 


(15) 


where  rn1  and  ra2  are  random  numbers  in  the  inter¬ 
val  (0,  1)  (uniform  deviates).  The  convolution  is  per¬ 
formed  for  a  suitably  large  number  of  anew  values, 
with  noise  added  to  each  convolved  55Fe  kernel. 
(Noise  values  are  randomly  selected  from  a  Gaussian 
distribution,  with  a  noise  standard  deviation  [<rnoise], 

ffnoise  =  1620G’  (16) 

where  G  is  the  measured  gain,  eq  is  the  expected  back¬ 
ground  noise  (in  e~),  and  1620  is  the  number  of  charge 
pairs  generated  by  a  Ka  photon  absorbed  in  silicon.) 
The  55Fe  kernels  are  summed,  and  the  final  set  of  pix¬ 
els  is  normalized  so  that  the  sum  of  charge  across  all 
pixels  in  the  final  kernel  =  1.  The  two-dimensional 
convolution  of  the  charge  diffusion  and  IPC  kernels 
is  performed  using  the  Python  routine  convolve2d.py. 
Figure  15  shows  the  results  of  the  55Fe  model. 

A  comparison  of  this  55  Fe  model  kernel  with  the 
measured  55  Fe  data  in  Fig.  14  indicates  that  the  com¬ 
bination  our  SPR  model  and  charge  diffusion  model 
predicts  the  actual  55  Fe  kernel  very  well  at  40  V.  The 
central  values  for  the  model  kernel  are  within  the  er¬ 
rors  of  the  actual  55 Fe  data,  as  both  tables  above 
show.  The  SPR  and  charge  diffusion  model  also  pre¬ 
dict  the  values  in  the  surrounding  pixels  very  well. 
We  note  that  some  of  the  pixels  around  the  central 
pixel  do  not  quite  match  the  model  to  within  the  er¬ 
rors.  This  is  likely  because  the  signal  in  the  sur¬ 
rounding  pixels  is  much  less  than  in  the  central 
pixel.  When  these  pixels  are  normalized  to  the  total 
in  the  kernel,  the  effect  of  noise  in  the  normalization 
is  much  higher  for  these  pixels  than  in  the  central 
pixel.  Thus,  in  addition  to  the  statistical  errors  across 


7  x  7  x  2419  Stamp  Stack 


OOOOOOL 

toQOOOQL 

HUBBUB. 


HHEEEHEBL 

Innmnri 

Innnnnnr 

n  n  n  n  n  n  — l 

_//  ll  ll  U  1 1 — 

n  n  n  n  n  x 


Fit  to  Pixel  Values  Across  2419  Stamps 


Fit  to  Value  Bins 
3  Sigma  Limit 


100  150 

Counts  in  Pixel 


Fig.  14.  (Color  online)  Kernel  and  errors,  with  pixel  values  resulting  from  55 Fe  tests  of  H4RG-10  B0. 
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Fig.  15.  (Color  online)  Predicted  kernel,  with  pixel  values  result¬ 
ing  from  the  55 Fe  model. 


all  hits,  there  is  error  in  the  normalization  due  to  the 
noise  of  the  system.  A  system  with  lower  noise  than 
the  camera  presented  here  will  likely  find  model  va¬ 
lues  even  closer  to  the  observed  values.  We  explore 
the  predictive  power  of  our  model  in  more  detail  in 
the  next  section. 

F.  55Fe  Modeling  across  Voltages 

55  Fe  measurements  were  taken  at  substrate  voltages 
ranging  from  10  to  40  V.  This  section  describes  the 
crosstalk  results  for  all  voltage  settings  and  com¬ 
pares  these  results  with  the  model  predictions. 

When  extracting  the  kernel  for  the  10,  20,  and  30  V 
measurements,  it  is  critical  to  take  into  consideration 
broader  spreading  induced  by  charge  diffusion.  For 
10,  20,  and  30  V  settings,  we  found  that  15  x  15, 
13  x  13,  and  11  x  11  square  pixel  kernels,  respec¬ 
tively,  were  needed  to  capture  all  of  the  charge  from 
the  55Fe  hit.  Because  the  gain  of  the  system  should 
not  change  with  substrate  voltage,  the  gain  from  the 
40  V  measurement  was  used  for  the  lower  voltage 
settings,  since  here  the  Ka  distribution  is  much 
tighter. 


Table  3.  55 Fe  Model  and  Data  Results  for  Substrate  Voltages  1 0-40  V 


Model  and  data  comparison  for  10-40  V  55  Fe 

Vsub  (V) 

Kernel  size 

Model 

Central  pixel 

Data  Std  dev 

10 

15  x  15 

32.24% 

29.16% 

2.61% 

20 

13  x  13 

47.63% 

46.85% 

2.79% 

30 

11  x  11 

55.57% 

54.65% 

1.65% 

40 

9x9 

60.51% 

60.43% 

1.58% 

X-right  pixel 

10 

15  x  15 

9.52% 

8.88% 

0.57% 

20 

13  x  13 

7.79% 

7.65% 

0.38% 

30 

11  x  11 

6.48% 

6.63% 

0.27% 

40 

9x9 

5.72% 

5.92% 

0.17% 

X-left  pixel 

10 

15  x  15 

9.50% 

8.88% 

0.25% 

20 

13  x  13 

7.77% 

7.67% 

0.41% 

30 

11  x  11 

6.48% 

6.57% 

0.36% 

40 

9x9 

5.71% 

5.94% 

0.16% 

Table  3  shows  the  central  and  neighboring  pixels 
resulting  from  the  55 Fe  measurements  and  those  pre¬ 
dicted  by  the  model. 

Figure  16  contains  plots  of  the  central  pixel  value 
and  nearest  neighbor  value,  both  for  the  measure¬ 
ments  and  for  the  model.  While  the  model  comes 
close  to  the  results  we  have  for  the  10  V  data,  it  does 
not  predict  the  data  to  within  the  errors.  At  10  V,  the 
distributed  charge  inhabits  a  very  large  kernel:  15  x 
15  pixels,  for  which  the  charge  is  very  close  to  the 
noise  floor  of  the  system.  We  therefore  conclude  that 
the  10  V  data  presented  here  is  not  high  fidelity. 

5.  Summary 

While  the  design  of  source-follower  CMOS-hybrid  ar¬ 
rays  eliminates  some  noise  sources  such  as  charge 
transfer  inefficiency  (CTI),  read  noise,  and  dark  cur¬ 
rent,  it  also  results  in  an  increase  in  IPC,  in  contrast 
to  other  detectors.  IPC  can  be  problematic  for  some 
astronomical  applications  such  as  astrometry  and 
photometry,  which  relies  on  high  S/N  to  obtain  accu¬ 
rate  and  scientifically  useful  results. 
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Fig.  16.  (Color  online)  Comparison  of  central  and  nearest  neighbor  (right  and  left)  pixel  values:  predicted  by  the  model,  and  taken  from 
the  final  55 Fe  kernel. 
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In  this  paper,  we  have  used  the  SPR  method  to 
model  IPC  in  TIS  H4RG-10  arrays.  We  combined  a 
charge  diffusion  model  with  this  IPC  model  to  under¬ 
stand  the  global  effect  of  crosstalk  on  these  arrays. 
Our  primary  findings  are  as  follows: 

•  The  IPC  in  the  H4RG-10  BO  detectors  extends 
to  9  x  9  pixels  based  on  our  SPR  method,  requiring  a 
extension  to  Moore’s  IPC  model  to  capture  all  of  the 
charge  spread. 

•  The  SPR  method  shows  a  readout  bias  that  is 
column-dependent.  This  effect  introduces  an  inher¬ 
ent  systematic  bias  in  the  IPC  pattern.  This  bias 
must  be  included  in  the  55  Fe  reduction  strategy  in 
order  to  extract  all  good  hits  from  a  55 Fe  image. 

•  We  have  developed  a  model  for  combining  IPC 
and  charge  diffusion  effects  as  a  function  of  substrate 
voltage.  Our  model  results  closely  match  our  empiri¬ 
cal  results. 

•  Our  results  suggest  that  the  SPR  method  alone 
combined  with  the  model  of  charge  diffusion  de¬ 
scribed  here  can  provide  enough  information  to  simu¬ 
late  the  dominant  properties  of  H4RG-10  BO  crosstalk 
to  high  accuracy  without  need  for  a  radioactive  source. 
In  addition,  the  SPR  data  with  this  charge  diffusion 
model  can  be  used  to  very  accurately  predict  astro- 
metric  and  photometric  performance  [4] . 

•  We  find  that  the  H4RG-10  BO  generation  of  de¬ 
tectors  shows  crosstalk  of  -10%,  and  a  signal  loss  of 
—40%  in  the  central  pixel  at  40  V.  A  risk  reduction  ef¬ 
fort  is  currently  underway  to  reduce  this  crosstalk  by 
half  in  a  new  generation  (Bl)  of  H4RG-10  detectors. 

The  model  based  on  SPR  data  presented  here  re¬ 
produces  the  55  Fe  data  to  within  the  statistical  errors 
and  noise  properties  of  our  system.  In  a  companion 
paper,  we  use  the  models  to  explore  the  effects  cross¬ 
talk  at  various  levels  (5-15%)  has  on  astrometric  and 
photometric  applications  [4]. 
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